Dielectric characterization and modelling of aqueous solutions involving sodium chloride and sucrose and application to the design of a bi-parameter RF-sensor

This paper reports on dielectric properties of ternary mixtures involving sodium chloride (NaCl) and sucrose (C12H22O11) dissolved into water (H2O). Broadband electromagnetic characterizations of such mixtures at various concentrations were performed, evidencing a dual behavior made of conductive effects at low frequencies and dipolar relaxation at microwave frequencies. Conductive and dielectric properties resulting from these both effects were integrated into predictive models for variations of Cole–Cole model parameters. Based upon this modelling, an innovative microwave-based sensor able to retrieve concentrations of both sodium chloride and sucrose in ternary aqueous solutions was introduced, designed, realized and assessed. The proposed sensor shows an error lower than 5.5% for concentration ranges of 0 to 154 mmol/L for sodium chloride and 0 to 877 mmol/L for sucrose.

Accordingly, techniques using low-volume samples are arousing attention, which lead to a lot of studies focusing on planar resonant techniques because of their ability to achieve measurements with relatively small samples 2,3,28 . In addition to the limited requested volume of liquid (typically few microliters), these planar resonant techniques draw attention thanks to their compactness, low cost and ease of manufacture. Furthermore, resonant RF-sensors usually exhibit greater sensitivities than broadband sensors. Several topologies can be found for resonant planar RF-sensors, from simple stub or LC resonator 29 to more complex designs such as Split Ring Resonators (SRR) [30][31][32] for example, which are quite convenient in sensing applications for combining the accuracy of the large structures while being miniaturized and highly sensitive. Although the sensing method of such sensors is usually based on the modification of the resonant frequency when in the presence of the sample to be characterized 3,15 , recent papers have shown that higher sensitivities can be reached thanks to the use of other indicators, such as the phase shift, magnitude of Scattering-Parameters or even Q-factors [33][34][35][36][37] .
Majority of those papers and the studies found in the literature are focused on the design of sensors able to track concentration of one given compound dissolved into a solution or a complex media. Some of them investigate or give clues about the cross-sensitivities of the tracked parameter with respect to other parameters such as temperature for instance. However, to our knowledge, very few papers report on the design of RF sensors dedicated to the tracking of several separate parameters with one single device in spite of the major interest it can provide for analyzing complex media such as those found in many applications.
A recently-published study focusing on sodium chloride and glucose sensing has been proposed 38 , it proposes to track concentrations of ternary mixtures through the visualization of microwave near-field distribution achieved thanks to thermo-elastic optical indication. As this work relies on a hybrid microwave/optical technique, it requires a quite complex set up and suffers from the optical methods classical drawbacks.
In the present paper, we propose to investigate the design of a single RF-based sensor purposed to determine the respective concentrations of two distinct compounds dissolved into an aqueous solution. Novelty of this paper is to give a proof of feasibility of such a dual-parameter sensor, its sensing abilities being exclusively based on electrical measurements. The proposed compounds to be tracked are sodium chloride and sucrose, which, as a first approximation, modify ionic behavior and dipolar relaxation, respectively. These intrinsic differences, together with the fact that these substances are widely used and found in plenty of applications have driven our choices. Indeed, sodium chloride and sucrose play a fundamental role in many chemical processes in a wide variety of domains.
Sodium chloride is an ionic compound, its molecular formula is NaCl. In its edible form of table salt, it is commonly used as a condiment and food preservative, it is largely found in nature as it is mainly responsible for the salinity of the seas. Large quantities of sodium chloride are used in many industrial processes, and it is a major source of sodium and chlorine compounds used as feedstocks for further chemical syntheses.
Sucrose is commonly known as table sugar, and is obtained from sugar cane or sugar beets, therefore, it is found in abundance in nature. Sucrose is a disaccharide formed from the monosaccharide glucose associated with a monosaccharide fructose. The molecular formula of sucrose is C 12 H 22 O 11 39,40 . Accurate detection and quantification of sucrose and sodium chloride concentrations in water with high sensitivity, selectivity and accuracy is thus of great interest. Notwithstanding the primary aim of this work was to demonstrate the feasibility of an RF-based sensor to track both sodium chloride and sucrose in aqueous solutions, concentration ranges consistent with those found in the food and beverage industry [41][42][43] were targeted (0-9 g/L, which stands for 0-154 mmol/L for sodium chloride, and 0-300 g/L, which stands for 0-877 mmol/L for sucrose) to highlight the applicative interest of this study.
Although the Electro-Magnetic (EM) properties of sodium chloride and sucrose are different, with ionic conductivity on one hand and exhibiting a dipolar relaxation on the other hand, these species interact when they are mixed together into a solution, especially through the modification of the viscosity of the solution. It is thus paramount to study these physical interactions and their impact over their respective dielectric signatures. To reach this objective, direct electromagnetic characterizations were performed. They allowed the obtention of a thoroughful EM models for ternary mixes involving water, sodium chloride and sucrose.
EM Characterization, analysis and model of binary and ternary mixes are presented in the next section of this paper.
Then, based on these EM models, and their specific features, a sensor was designed and realized, which is presented in "Sensor design" section. "Experimental assessment of the sensor" section is devoted to the measurements and assessment of the fabricated sensor together with discussion of the obtained results. Finally, conclusions and future scopes are proposed in the last section.
Dielectric characterization and modelling of aqueous solutions. The polarization undergone by a material when an electrical field is applied is referred as dielectric permittivity; it is generally given as a complex value: where ω is the angular frequency and j the imaginary unit.
The real part ε ′ of the dielectric permittivity is a measure for the polarizability of a material and the imaginary part ε ′′ , is a measure for dielectric losses 25,44,45 . In aqueous solutions, the polarizability originates from rotation of the dipole moment of water molecule and the large value of ε' comes from cooperative motions due to intermolecular hydrogen bonds 46,47 . The dielectric losses ε" mainly arise from the dipole relaxation, which occurs at microwave frequencies. Moreover, for aqueous solutions with dissolved salts, additional conductive losses arising from ion mobility occur at low frequencies. The most commonly used model to describe the complex where ε ∞ and ε stat are the limits of the relative dielectric permittivity at low and high frequencies, respectively, ω is the angular frequency, τ is the relaxation time, σ s is the static ionic conductivity, ε 0 is the free space permittivity (8.854 × 10 −12 F/m) and α is a dimensionless exponent parameter taking into account dispersion in relaxation time. For De-Ionized (DI) water at room temperature (23 °C) the commonly accepted values of Cole-Cole parameters 49 are: ε ∞ = 5.1, ε stat = 79.2, τ = 8.72 ps, α = 0 and σ S = 0 mS/cm. Most aqueous solutions with dissolved solutes can thus be modelled using Eq. (2), whose parameters depend on solute concentrations and other quantities such as temperature or pressure.In this work, we have chosen to study two common solutes that influence the behavior of the permittivity in different ways. Firstly, sodium chloride was used, it adds dissolved ions and mainly changes the conductive behavior of the solution. Secondly, sucrose was used, it increases the viscosity and mainly modifies the dipolar behavior of the solution. In next parts, we propose to extract the concentration dependency of parameters of Eq. (2) at room temperature and pressure; firstly for sodium chloride solutions, then for sucrose solutions and finally for mixtures of both sodium chloride and sucrose.
To do this, broadband permittivity measurements were performed on aqueous solutions using the openended coaxial probe method (Keysight ® N1501A-104 probe and Agilent ® E8364 PNA-L Network Analyzer) 27 in the [100 MHz-50 GHz] frequency range at room temperature (23 °C). A Matlab ® curve fitting code was developed to extract values of Cole-Cole model parameters that best fit the measured complex permittivity. Taking advantage that the variations of the imaginary part at low frequency are important but mostly depend on the conductive effects we chose a two-step fitting procedure to reduce optimization uncertainties and improve the accuracy of optimized parameters ε ∞ , ε s , τ, α and σ s . First, the conductivity σ s is determined from fitting the imaginary part of the permittivity at low frequency (Fig. 1a). Then the conductive behavior is subtracted from the measurement according to Eq. (2) to keep only the dipolar behavior. Finally, the model parameters ε ∞ , ε s , τ and α are determined from fitting both imaginary and real parts of the dipolar permittivity (Fig. 1b). In addition, measurements of static ionic conductivity were performed using a SevenExcellence™ EC meter from Metter Toledo to validate the method. Agreement between static measured conductivity values and those determined from dynamic permittivity is better than 2%. Figure 1 shows an example of the two-step fitting procedure to extract the model parameters from permittivity measurement of a salt water solution.
Note that in Eq. (2), ε ∞ is a mathematical parameter that represents the dielectric permittivity at high frequencies but does not correspond to the actual physical value of permittivity in the infrared frequency range 49,50 . The refractive index n = √ ε ∞ calculated from the infinite permittivity of the Cole-Cole model (ε ∞ = 5.1) is different from the 1.33 measured by optical techniques. Other relaxation and resonance processes appear at high frequencies that are not taken into account in the Cole-Cole model 49 . The infinite permittivity in the Cole-Cole  Although the value of ε ∞ is sometimes fixed to the water one in the literature 51 , we decided to use it as an optimization parameter to best fit the experimental curves in our frequency range of interest without relating its variation to physical phenomena.
Sodium chloride aqueous solution. Due to its large polarity, water is an excellent solvent that is able to dissociate salts by separating the cations and anions and forming new interactions between the water and ions. The dielectric properties of sodium chloride aqueous solutions have been intensively studied in the literature [51][52][53] . The sodium chloride dissolves to sodium cation and chloride anion, thus providing an ionic character and helping to make the solution conductive.
In the presence of an electric field, ions move according to their mobility, which strongly depends on the viscosity 51 of the solution. Indeed, ionic conductivity is proportional to the concentration of sodium chloride in the solution and inversely proportional to the viscosity [54][55][56] . It is worthwhile to note that the temperature influences the viscosity, this is highlighted in various research works 28,57,58 .
In addition to ions that generate a conductive effect, the sodium chloride aqueous solutions also exhibit a dipolar relaxation alike the DI water one. This dipolar relaxation being related to the viscosity of the solution, it is also influenced by temperature and solute concentrations. To evidence these effects on EM properties, permittivity characterization of sodium chloride aqueous solutions was performed.
Solutions were made by diluting sodium chloride from Merck ® (6404) into DI water. Sodium chloride was weighted thanks to a precision balance from Kern Tab ® . Obtained concentrations of sodium chloride in the solution were respectively 51 mmol/L, 103 mmol/L and 154 mmol/L. The measured complex permittivity of sodium chloride aqueous solutions is shown in Fig. 2.
For non-zero sodium chloride concentrations, the imaginary part of the permittivity exhibits two regions of losses. First at low frequencies (up to 1 GHz approximately), we observe a decrease in the imaginary part inversely proportional to the frequency and proportional to the sodium chloride concentration corresponding to the ionic conductive losses. Then at high frequencies, we observe the losses of dipolar relaxation of water molecules with a maximum of the imaginary part occurring around 18 GHz whatever the concentration.
The real part of the permittivity slightly decreases as sodium chloride concentration increases. For high concentrations, measurement uncertainties appear at low frequency due to the electrode polarization effects 59 . So the first frequency points of measured permittivity were ignored for the curve fitting of high concentrations.
From these measurements, we observe, as expected, that sodium chloride clearly modifies the ionic behavior of the permittivity and hardly modifies the dipolar behavior. The extracted parameters of the model are given in Table 1 for each concentration.
From this table, one can note that the static permittivity ε stat is slowly decreasing as the sodium chloride concentration increases. This is the consequence of, first, the lower concentrations of water when sodium chloride concentration increases and, second, of the ion-solvent interaction which leads to non-rotational bonding of water molecule in the vicinity of ions.  www.nature.com/scientificreports/ The relaxation time τ slightly decreases as sodium chloride concentration increases. This behavior is quite unexpected because increasing the solute concentration generally leads to increase the viscosity and thus to decrease the relaxation time. However, in sodium chloride solutions, viscosity increases very slightly with concentration 60 . Moreover, two antagonist contributions are present, a relaxation time of water surrounding cation (τ+) higher than those of DI water and a lower relaxation time of water surrounding anion (τ−) 61 . The summation leads to a slightly decreasing mean relaxation time for sodium chloride concentrations below 1 mol/L 51 .
The ionic conductivity σ s considerably increases as the sodium chloride concentration increases. As the viscosity of sodium chloride solutions shows very little changes, the conductivity behaviour is mainly driven by the amount of ions added to the DI water.
The parameter α slightly increases as the sodium chloride concentration increases. The presence of solutes in solution leads to a broadening of the relaxation losses that is well described by dispersion in relaxation time.
The parameter ε ∞ decreases slightly as the concentration increases but since it acts as a mathematical optimization parameter, no physical phenomena was associated to its variation.
The parameters extracted from our measurements on sodium chloride aqueous solutions (Table 1) were compared to those found in the literature 51 . The values we obtained and their variations with solute concentration are in agreement with previously-published models.
To summarize, sodium chloride diluted in water mainly provides a conductive character due to dissolved ions, and only slightly affects the dipolar character due to water molecules. Note that the concentration of sodium chloride can be determined by knowing only one parameter such as the static ionic conductivity of the solution. Now let us investigate the dilution of a non-ionic solute in water, and its effects on the Cole-Cole model parameters.
Sucrose/water mixtures. Sucrose molecule was chosen because it greatly affects the viscosity of aqueous solutions and does not exhibit ionic character. When diluted into water, the dipolar character of the solution mainly comes from the water molecules, but it is greatly affected by the viscosity modifications brought by sucrose. There is a large volume of published studies describing the viscosity properties of sucrose aqueous solutions 58,62 . However, there are relatively few studies in the area of permittivity modelling of sucrose solutions. A work 63 studied the permittivity dependency on solute concentration, but this study was limited to static case only. Another published work 64 proposes the study of the complex permittivity of carbohydrate aqueous solutions, but only for one concentration (33% w/w = 1.1 mol/L) at a low temperature of 5 °C. So, to evidence effects of sucrose on EM properties, permittivity characterization of sucrose aqueous solutions was performed.
Solutions were made from pure sucrose dissolved into DI water. Sucrose was weighted thanks to a Precisa XB1200C precision balance. Obtained concentrations of sucrose in the solutions were respectively 293 mmol/L, 585 mmol/L and 879 mmol/L. Figure 3 shows the permittivity measurement of sucrose aqueous solutions at various concentrations over a frequency range starting at 100 MHz and up to 50 GHz.
As expected, no conductivity effects on these measurements of permittivity can be seen. The imaginary part of the permittivity exhibits only one region of losses at high frequencies corresponding to the relaxation process. However, the shift towards low frequencies of the maximum of the imaginary part shows the effect of viscosity brought by sucrose. Moreover, the maximum values of both real and imaginary parts of the permittivity significantly decrease as sucrose concentration increases. Table 2 summarizes the values of model parameters extracted from measurements at each concentration of sucrose.
The static permittivity ε stat significantly decreases as the sucrose concentration increases. As for sodium chloride, this is a consequence of the lower concentrations of water molecules when sucrose concentration increases. Moreover, the sucrose molecule is large and its first hydration shell contains over 37 water molecules 65 leading to much more non-rotational bonds than for sodium chloride. That is why the decrease of the static permittivity The relaxation time τ is greatly impacted by sucrose concentration because it strongly affects viscosity 58 , reducing water molecules mobility and thus their ability to relax quickly when the dynamic electric field is applied. Thus, the relaxation time significantly increases as the sucrose concentration increases.
The static conductivity σ s is null as sucrose does not introduce any ionic contribution in aqueous solutions.
The parameter α increases significantly as the sucrose concentration increases. The widening of the imaginary part of the permittivity spectra is greater than for sodium chloride solutions because the viscosity changes brought by sucrose are greater.
As for sodium chloride solution, the parameter ε ∞ decreases as sucrose concentrations increase.
We have shown that, as expected, sucrose significantly modifies the dipolar behavior and does not modify the conductive behavior of the permittivity. Note that the concentration of a sucrose aqueous solution can be determined by knowing only one parameter such as the static permittivity of the solution.
Now that effects over EM characteristics of sodium chloride aqueous solutions on one hand and sucrose aqueous solutions on the other hand, have been described, we propose in the next section to investigate the behavior of ternary solutions involving both sodium chloride and sucrose diluted into water.
Ternary mixtures involving water, sodium chloride and sucrose. In ternary solutions made from sodium chloride, sucrose, and water, several physical phenomena occur and interactions between them appear. Especially, the viscosity brought by sucrose significantly affects the ionic conductivity brought by sodium chloride. Although different studies exist in the literature regarding the viscosity or conductivity of such ternary solutions, none deals with the change of permittivity according to interactions between sucrose and sodium chloride in aqueous solutions.
To identify and quantify these interactions, broadband (100 MHz-50 GHz) permittivity measurements of ternary mixtures made up from sodium chloride, sucrose and water were performed at room temperature (23 °C). Solutions were prepared using the same compounds, the same procedure and the same equipment as in the previous sections. Targeted solutes concentrations are close to 0 mmol/L, 50 mmol/L, 100 mmol/L and 150 mmol/L for the sodium chloride and 0 mmol/L, 290 mmol/L, 580 mmol/L and 870 mmol/L for the sucrose concentrations.
To better identify the solutions, one letter is attributed to each sodium chloride concentration whereas one number is associated to each sucrose concentration. Thereby, A0 stands for DI-water, samples labelled B0, C0 and D0 are binary mixtures made from sodium chloride and water, as studied in "Sodium chloride aqueous solution" section and samples labelled A1, A2, A3 are binary mixtures made from sucrose and water, alike those presented in "Sucrose/water mixtures" section. All other labels with letter different than A and number different than 0 are ternary mixtures. Figure 4 shows the variation of measured complex permittivity of ternary mixtures over the frequency range 100 MHz-50 GHz.
As in the previous sections, we use the curve fitting method to extract the model parameters from these measurements. Table 3 shows the resulting parameters.
We observe that the static permittivity ε stat decreases when either the sodium chloride or the sucrose concentrations increase, because of the same reasons as for binary mixtures. In ternary mixtures, the decrease of ε stat is mainly due to the concentration of sucrose but is emphasized by the sodium chloride concentration.
As for binary mixtures, the relaxation time τ of ternary mixtures is strongly impacted by sucrose due to its strong impact on viscosity. Thus, the relaxation time increases significantly as the sucrose concentration increases but is less influenced by sodium chloride. For a fixed sucrose concentration, the relaxation time decreases slightly as the sodium chloride concentration increases for the same reasons as those given in "Sodium chloride aqueous solution" section.
Viscosity effects brought by sucrose are greatly noticeable on the ionic conductivity. For example, by considering a sodium chloride concentration set at 154 mmol/L, the conductivity σ S decreases from 15.99 to 8.12 mS/ cm as the sucrose concentration increases from 0 to 877 mmol/L respectively. The conductivity σ s of ternary solutions is mainly related to the amount of solvated ions: the greater the concentration of sodium chloride the greater the conductivity. However, viscosity reduces the ion mobility: the higher the sucrose concentration, the lower the conductivity.
As in binary mixture, the greater the concentration of sucrose in ternary mixture, the greater the α parameter. In addition, α slightly increases as the sodium chloride concentration increases. Thus, the higher the solute concentration in ternary mixture, the wider the loss curves.
Concerning parameter ε ∞ , the greater the concentration of sucrose, the lower the parameter. In addition, ε ∞ , decreases slightly as the sodium chloride concentration increases. www.nature.com/scientificreports/ In general, solutes interactions take place in ternary solutions and modify the variation laws of the parameters of binary mixtures. Variations of parameters ε s , τ , ε ∞ , α and σ s were plotted in Fig. 5 with respect to the concentrations of sodium chloride for the different concentrations of sucrose.
All parameters show quite linear variations depending on the sodium chloride concentration, slopes and intercepts of which depend on the sucrose concentration.
Now that the interaction effects occurring in ternary mixtures have been identified and model parameters have been extracted from measurement results, the next section will provide predictive numerical models for the variations of the model parameters regarding sodium chloride and sucrose concentrations.
Variational laws of the Cole-Cole model parameters for ternary mixtures. Polynomial equations were computed from the data reported in Table 3, relating each of the Cole-Cole model parameter (static permittivity, relaxation time, ionic conductivity, distribution parameter and infinite permittivity) to both the concentrations of   www.nature.com/scientificreports/ sucrose and sodium chloride. The following equations are provided under the assumption that all the concentrations are given in mol/L, the concentrations are within the 0 to 0.154 mol/L range for sodium chloride and 0 to 0.877 mol/L for sucrose, temperature was stable and set to 23 °C. In these conditions, static dielectric permittivity ε stat can be expressed as: where ε stat DIW is the static dielectric permittivity of DI-water. The absolute error between modelled and measured static permittivity is less than 0.2. Since static ionic conductivity strongly depends on the concentration of sodium chloride, we have chosen to model it precisely using a third-degree polynomial. So σ S can be expressed as: where τ is expressed in ps. The absolute error between modelled and measured relaxation time is less than 0.13 ps. Dispersion parameter α expression is: The absolute error between modelled and measured value of parameter α is less than 0.008. Optimization parameter ε ∞ can be expressed as: The absolute error between modelled and measured value of parameter ε ∞ is less than 0.2. Now that the model parameters have been extracted and their variations modelled as a function of sucrose and sodium chloride concentrations, the EM response of the mixtures can be predicted whatever the solutes concentrations in the given ranges. Furthermore, by solving the inverse problem, it is possible to determine sodium chloride and sucrose concentrations in a ternary mixture as long as at least two parameters of the model are known: the static permittivity and the ionic conductivity, for example. This is the starting point to develop microwave sensors dedicated to the monitoring of concentrations in ternary mixtures, which is detailed in the next section.
Sensor design. Sensing principle. As demonstrated before, extraction of concentrations can be made from Eqs. (3) and (4) as long as conductivity and real part of the static dielectric permittivity are known. However, from a single low-frequency measurement, separating effects of conductivity from those related to static dielectric permittivity is not possible without additional information: it is compulsory to obtain complementary data. In this paper, we propose to obtain this complementary information from microwave-range measurement. Indeed, as shown in the previous section, solute concentrations of ternary mixtures can be identified from microwave frequencies dielectric properties of the solutions. To achieve high sensitivities, choice was made to rely on a RF planar resonant sensor topology. In addition to great sensitivities that could be achieved, resonant sensors usually require relatively simple driving electronics and systems to perform measurements, enabling their easy integration and featuring nice portability potential.
To limit the handling of samples and to keep the number of items reduced to the minimum, the use of a single sensor employing a single sample to perform both RF and low-frequency measurements is highly advised. The idea was thus to develop a unique sensor designed to extract low-frequency conductivity and complementary information from RF-range to extract simultaneously sodium chloride and sucrose concentrations out of ternary mixtures.
To achieve conductivity measurements at low frequencies, a classical way is to design a serial capacitance between two ports and to place the Material Under Test (MUT) upon it, in contact with the electrodes. MUT affects the serial capacitance, which is thus short-circuited due to the presence and displacement of ions if the MUT is conductive. Measurement of the resulting resistivity enables extraction of the conductivity. Figure 6 shows the principle of such a conductivity sensor operating at low frequencies.
In this Figure, C stands for the serial capacitance whereas R σ , serial resistance, is modelling the conductive effect brought by the ions existing in the solution, R σ is proportional to angular frequency ω divided by the static conductivity σ s . Thus, R σ is very small at low frequencies (short-circuit behaviour) and is rapidly increasing as the frequency increases (almost open-circuit behavior).
Using such a serial-capacitance sensing principle for low frequencies enables the achievement of conductivity measurements, however, it also limits the possibilities of design for the microwave sensors.
Indeed, it means that (a) the only capacitance seen at low frequencies between input and output ports should be the one under the MUT, and (b) the design should be made so that it has great interactions at RF frequencies with the MUT located over the capacitance.
Therefore, all the topologies relying on capacitively-coupled accesses or those showing direct DC input-tooutput paths should be discarded.   www.nature.com/scientificreports/ Considering these design constraints, for the RF design, sensor topology presented in Ref. 66 was chosen as a starting point. This sensor is composed of two facing Split Ring Resonators (SRR), also called Open-Loop Resonators (OLR), with the sample placed in the electric coupling slot area. This topology was chosen because it exhibits high sensitivities to dielectric permittivity changes such as those of glucose, and consequently those of sucrose. Furthermore, this topology is compatible with the low frequency requirements described hereinbefore.
Sensing principle of such a topology can be modelled with the lumped-element circuit presented in Fig. 7.
In this model, G 1 , C 1 , L 1 and G 2 , C 2 , L 2 are the lumped elements modelling the first and second resonator, respectively. C stands for the coupling between the two resonators and finally G 1Sample , C 1Sample , G 2Sample , C 2Sample , R Sample and C Sample refer to the dielectric contributions of the sample (MUT) to each constitutive element of the sensor model. R σ models the conductive effect of the sample.
The sensor is thus built from two resonators coupled to each other, which are fed by input-output lines directly connected to the resonators. The sample is placed in the coupling area, at the open-ended tips of the split ring resonators (SRR). This area is chosen as it is the location where the electric field is maximum, thus the sample greatly affects the inter-resonator coupling (contributions of R Sample and C Sample ). On the other hand, it also affects the inner capacitance and resistance of each resonator as the sample is close to the open-ended tips of the resonators. This latter contribution of the sample over each resonator is modelled thanks to G 1Sample and C 1Sample for the first resonator and G 2Sample and C 2Sample for the second resonator.
One should note that if the physical implementation is strictly symmetrical yielding identical resonators, values with indices 1 and 2 are equal ( G 1 = G 2 , C 1 = C 2 , L 1 = L 2 , G 1Sample = G 2Sample and C 1Sample = C 2Sample ).
It is also worthwhile to note that for very low frequencies, resonators are not resonating, and this model can be simplified, to a model alike the one presented in Fig. 6. Reversely, model in Fig. 7 is relevant for a microwave sensor under the assumption that R σ is high enough (open-circuited behaviour) to ensure a good sensitivity with respect to R Sample and C Sample . This is related to the choice of the working frequency and the design geometry, which will be presented in the next section. Once the frequency range is chosen, variations in the dielectric properties of the samples ( ε ′ and ε ′′ ), which are related to sodium chloride and sucrose concentrations, will modify G 1Sample , C 1Sample , G 2Sample , C 2Sample , R Sample and C Sample , and consequently the associated electrical response of the sensor, which is observed through its scattering parameters.
Sensor design. Once the sensing topology is chosen, the next step toward the design of the sensor is to determine the targeted frequency range of operation for the microwave part. Different effects have to be considered to do so.
On one hand, microwave operating frequency should be high enough to evacuate the conductive effects of sodium chloride as shown in Fig. 2, which means that frequencies lower than 5 GHz should be discarded.
On the other hand, high frequencies (> 10 GHz) are not suitable because the real part of the permittivity for all water-based solutions is decreasing toward high frequencies. Using those high frequencies would mean low values for the real part of the permittivity (see Fig. 4), which would lead to less concentration of the electric field distribution within the sample, which is not favorable for sensing purposes. Additionally, the imaginary part of the dielectric permittivity is reaching a maximum around 12 GHz, which means increased losses at these frequencies.
To identify the best frequency range for sensing purposes, maximum variations of the real part of the dielectric permittivity obtained thanks to the wideband EM characterizations presented in "Ternary mixtures involving water, sodium chloride and sucrose" section were computed. Indeed, sensitivity of the sensor is mainly linked (although not exclusively) to the variations of C 1Sample , C 2Sample , and C Sample , which are mainly related to changes of the real part of the permittivity of solution.
In order to identify the best operating frequency, the difference between the permittivity real part of several solutions and those of DI water (DIW) was measured at each frequency for all solutions. The resulting difference ( �ε ′ f = ε DIW f − ε {Samples} f is plotted as a function of frequency (f) in Fig. 8 for three solutions.
This plot shows that maximum difference of the real part of the permittivity is reached around 8 GHz. www.nature.com/scientificreports/ Similarly, G 1Sample , G 2Sample , R Sample are mainly linked to the values of the imaginary part of the permittivity (Fig. 3), which should be chosen as low as possible to minimize losses, but whose variations might also be interesting for sensing purposes.
To trade-off between all these effects, an operating frequency around 7 GHz was chosen. At this frequency, the real part of the dielectric permittivity can be expressed as: where ε DIW (7GHz) is the real part of the dielectric permittivity of DI water at 7 GHz and ε 7GHz stands for the real part of the dielectric permittivity, at 7 GHz, of a mixture with a given sodium chloride concentration (C NaCl ) and a given sucrose concentration C C 12 H 22 O 11 in mol/L, all the other parameters and assumptions being the same as in "Ternary mixtures involving water, sodium chloride and sucrose" section. The absolute error between modelled and measured real part of permittivity at 7 GHz is less than 0.2. At this point, it is worthy to note the increase of sensitivity to sucrose concentration at 7 GHz with respect to that presented in Eq. (3) at DC, reinforcing the relevancy to choose RF range for the operating frequency.
To implement the sensor, a 0.8 mm-thick RT Duroid 5880 substrate with a relative dielectric permittivity of 2.2 and a loss tangent of 0.0009 was chosen.
Concerning the sample, it must be placed in the area where the electrical field is maximum, i.e. in the slot between the two resonators. Additionally, the liquid sample has to be in contact with the metal strips to allow for conductivity measurement at low frequencies. Consequently, a bottomless sample holder was made thanks to a 3D-printer. This annular-wall shaped piece was then stuck onto the designed circuit. Material used to achieve the sample-holder is RS-F2-GPWH-04 from FormLabs ® , once cured, its expected relative dielectric permittivity is 3.55 and loss tangent is 0.01. The glue to stick the sample holder and ensure a proper sealing is Loctite™ EA 3421from Henkel ® .
A simulated operating frequency around 7 GHz when loaded with samples was achieved with dimensions of L = 10 mm and l = 5.4 mm. The remaining parameters (w, width of the SRR lines, equal to 1.8 mm), (g, slot between the two resonators, equal to 1.2 mm) and (e, open-end gap of the SRR equal to 2 mm) did show remarkable influence upon the sensitivities (both to sodium chloride and sucrose concentrations), and had to be tuned in a recursive manner to obtain the best trade-off. Finally, input/output lines were designed to show a characteristic impedance of 50 ohms (T L width = 2.474 mm), their location being tuned to insure an appropriate coupling (Z = 5 mm).
The layout of the resulting sensor and its dimensions are shown in Fig. 9.
Size of the bottomless sample holder was tuned to R = 3.4 mm × r = 2.04 mm, which correspond to an estimated sample volume of 43.56 μL. The MUT container is located in the area where the electric field is maximum, and centered to ensure a symmetrical design.
Full wave simulations (run with HFSS™ software from ANSYS ® ) of the sensor with air and DI-water as MUTs are shown in Fig. 10. One should note that DI-water was modelled at this stage with the Cole-Cole model parameters extracted from the above-mentioned wide-band characterization of the A0 sample. Consequently, in the simulations, the DI-water, considered at room temperature, was modelled with a static permittivity ε stat of 79.2, an infinite permittivity ε ∞ of 4.6 and a relaxation time τ of 8.71 ps.
Adjustment of the sensor model. The sensor parts were realized and assembled; a photography of the realized sensor is shown in Fig. 11.
The sensor was then measured using a properly-calibrated Vector Network Analyzer (Rhode & Schwarz ® ZVA 67). Simulations and measurements of the sensor with an empty container (i.e. air as MUT) are shown in www.nature.com/scientificreports/ Fig. 12a together with the simulated and measured results when the container is filled with DI-water (A0 solution), which are presented in Fig. 12b. As it can be noted from Fig. 12, simulations of the initial model of the sensor are slightly different from measurements. It means that the proposed 3D-model of the sensor does not fit exactly the measurements of the realized sensor, the differences, mainly a frequency shift, can be attributed to inaccurate material description or technological dispersions. Concerning those latter ones, they can be related to machine tolerances or deviations, which could impact metal etching or container size. Technological dispersions may also be due to the manmade gluing of the container, entailing uncertainties upon the size of the final glue tie and the container position.
To consider all those dispersions, retrofitting electromagnetic simulations were run to refine the sensor model trying to make its simulations fit both the air and DI-water measurements. Sensitive design parameters were identified and tuned. The modified dimensions are g = 1.6 mm and size of the bottomless container R = 3.31 mm   www.nature.com/scientificreports/ × r = 2.06 mm. The glue tie model and dimensions were also slightly modified, relative changes on their parameters being lower than 5%. Concerning the material properties tuning, materials of the sensor only were tuned, MUT models (air or DI-water) were kept untouched during this retrofit step, and simulations yield a change on the relative dielectric permittivity of the substrate only, moving from 2.23 to 2.1.
Finally, the corrected model simulations are plotted in Fig. 13 together with the measured results.
Evaluation of sensor sensitivities from simulated data. To determine the sensing capability of the proposed design, full-wave simulations of the updated sensor were run in HFSS™ with various samples. These latter were modelled thanks to a Cole-Cole model with the parameters provided in "Dielectric characterization and modelling of aqueous solutions" section (Table 3). Broadband simulations were conducted but two frequency ranges were particularly scrutinized: the first being at low frequencies (between 10 MHz and 2 GHz) to track conductive effects while the second is located at higher frequencies (around the 7 GHz), following explanations given in "Dielectric characterization and modelling of aqueous solutions" section. Figure 14 shows the simulated transmission S-parameter for the sensor with various modelled samples. One can note the distinctive effects at low frequencies occurring when the modelled samples exhibit sodium chloride concentration changes (modelled samples A0, B0, C0, D0) whereas, as foreseen, samples modelling changes in sucrose concentration impact the transmitted S-parameter in the microwave range.
From these EM simulations, and aimed to assess sensor performances, several sensing magnitudes could be chosen, resulting in a trade-off between sensitivity, ease of measurement, data extraction and computation easiness, and error level related to measurement noise floor 67 .  www.nature.com/scientificreports/ Concerning low frequencies, since the objective is to extract the conductivity, sensing magnitudes could be either the magnitude of the transmitted S-parameter (S 21 ) or the reflected S-parameter (S 11 ). We chose the latter one (in natural value) because of its higher level and better sensitivity, and the frequency of 100 MHz was chosen to extract this sensing parameter.
In the microwave range, as detailed in Ref. 66 , several indicators could be chosen as sensing parameters, such as for instance, magnitude of the transmission S-parameter (S 21 ) or coupling coefficient K. The sensor being composed of two coupled resonators, its electrical response is alike a second order filter, and following the design choices detailed in "Sensor design" section, the coupling between the two resonators is strong, preventing the appearance of two distinct peaks in the transmitted scattering parameter. As the sample affects not only the mutual coupling (C Sample and R Sample , see Fig. 7) and consequently the bandwidth and insertion losses of the second order filter, but also intrinsic resonant frequency and Q-factor of each resonator through C 1Sample , C 2Sample , R 1Sample , and R 2Sample , we chose to define a sensing parameter, alike an unloaded Q-factor, but for a second-order filter. Definition of this dimensionless factor called Sensing Ratio (SR) is given hereafter: where BW is the bandwidth computed from 1 dB fall from S 21 maximum (instead of the commonly-used 3 dB fall, which cannot be computed for high sucrose concentrations, because of losses), f c is the central resonant frequency and S 21 f c is the magnitude (in natural value) of the transmission S-parameter at the central frequency. Table 4 summarizes values of those two sensing parameters ( |S 11 (100 MHz)| & SR) for simulations performed with each sample modelled with the parameters provided in Table 3.
From Table 4, simulated value of |S 11 (100 MHz)| for each sample was subtracted from that of the water (A0, with a |S 11 (100 MHz)| = 0.996 ) and plotted as a function of sodium chloride concentration. The resulting plot can be seen in Fig. 15a, where sodium chloride concentration is given in x-axis and expressed in mmol/L. At this point, one should note that this is fully consistent with the behavior related to the conductivity in the measured samples evidenced in Fig. 5c, demonstrating the relevancy of the design sensor on one hand and of the chosen sensing parameter on the other hand. Figure 15a shows the impact of sucrose concentration over the |S 11 (100 MHz)| sensitivity, which could also be considered as cross-sensitivity of sucrose concentration over sodium chloride concentration sensitivity of the sensor. Polynomial equations were computed from the data reported in Fig. 15. The following equations are provided under the assumption that all the concentrations are given in mol/L and the temperature is stable and set to 23 °C.
Since values of |S 11 (100 MHz)| strongly depends on the concentration of sodium chloride, we have chosen to model it precisely using a second-order polynomial. When standardized with respect to that of DI water, |S 11 (100 MHz)| of a given sample can be expressed as: (10) |S 11 (100 MHz)| = |S 11 (100 MHz)| DIwater − a i C 2 NaCl + b i C NaCl .   Fig. 15b as a function of sucrose concentrations. This plot shows the impact of sucrose concentration over the |S 11 (100 MHz)| sensitivity. Again, a linear approximation was made, which led to a regression coefficient of R 2 = 0.966, showing a good linearity. Expression of these coefficients are: The same approach was then performed concerning SR from Table 4, the computed SR value for each modelled sample is subtracted from that of DI-water (for A0, SR = 9.818). These SR standardized with respect to that of DI-water are then plotted as a function of sucrose concentrations for various sodium chloride concentrations in Fig. 16a. In this Figure with a regression coefficient better than 0.9987. Slopes of each line represent the simulated sensitivity of SR of the designed sensor to sucrose concentration variations. Figure 16 shows that the Sensing Ratio, which is alike a Q-factor computed at microwave frequencies from the S21 level, the central frequency and the bandwidth, is decreasing when sucrose concentration increases. This can be explained through the inverse behavior with respect to the relaxation time τ and the dispersion parameter α , which are both increasing when sucrose concentrations increase (Figs. 4, 5b,d).
Values of the slopes (coefficient c i ) and origins (coefficient d i ) are summarized in Fig. 16b as a function of sodium chloride concentration. This plot shows the impact of sodium chloride concentration over SR sensitivity, which represents the cross-sensitivity of sodium chloride concentration over sucrose concentration sensitivity for the designed sensor.
A linear approximation with a regression coefficient of R 2 = 0.995 was made. Expression of these coefficients are: With sodium chloride concentration (C NaCl ) and sucrose concentration C C 12 H 22 O 11 expressed in mol/L. Now that sensing magnitudes |S 11 (100 MHz)| and SR have been defined and their sensitivities with respect to sodium chloride and sucrose concentrations have been analyzed from simulated results, next section will focus on the retrieval of sample concentrations from measured electrical responses of the sensor, which relies on solving the inverse problem for Eqs. (10) and (13).   (Table 5). After each sample measurement, a careful double rinsing procedure based on ethanol and DI-water was applied. It was followed by a compressed-air drying step before moving on to the next sample measurement. All measurements were made at a room temperature of 23° C. Measurement results are presented in Fig. 17.
At this point, one should note that broadband measurements are shown for the exhaustivity sake, but retrieval of the concentrations are based on the measured S-parameters at 100 MHz on one hand, and on a frequency band located between 6 and 8 GHz on the other hand. Table 5 presents the values of the return loss level at 100 MHz ( |S 11 (100 MHz)| ) and the Sensing Ratio (SR, computed from Eq. (7) definition) directly extracted from the measurements shown in Fig. 17.
Then, from those measured values of |S 11 (100 MHz)| and SR on one hand, and Eqs. (10)-(15) on the other hand, solving the inverse problem yields retrieved concentrations for sucrose and sodium chloride. Those retrieved concentrations, extracted from measurements, are also reported in Table 5.
This table shows that it is possible to extract both sucrose and sodium chloride concentrations from the measurement of a single sensor. However, in order to further evaluate the error made while computing the concentrations and thus the overall sensor accuracy, those retrieved concentrations have to be compared to the real ones. Table 6 summarizes, for all samples, effective concentrations (labelled "Real Concentrations", which are the same data than those in Table 5) and retrieved concentrations (i.e. obtained by solving the inverse problem of Eqs. (10)-(15) upon measured |S 11 (100 MHz)| and SR). One should note that the inverse problem solving has led to unrealistic negative (yet very close to zero) values for the concentrations, which have consequently been set to zero. The relative differences for each compound concentration together with the quadratic error were also computed and reported in the same table.
Quadratic error is computed thanks to the formula: Analysis of this table shows that the extracted values for concentrations are consistent with the real ones, the relative error on each compound concentration being smaller than 9.8% for the sodium chloride and 9.5% for the sucrose. Average quadratic error over the set of measured samples is 11.1%.
These results show that the study proposed in this paper enables the detection of two distinct compounds in aqueous mixtures with a single sensor, yielding an overall error over species concentrations that does not exceed 11.1%. This validates the proposed models and the choices made on the scientific approach: modelling, sensor design and sensing parameters. This also substantiates the concept of an RF-based bi-parameter sensor, which has never been previously proposed, to the best of our knowledge, in the scientific literature. However, errors presented in Table 6 could be detrimental to industrial applications where high accuracies are seeked and consequently, these errors or inaccuracies should be analyzed to identify their causes.
As described hereinbefore, the proposed study involves numerous steps and modelling stages, each of them introducing potential errors.
To evaluate the accuracy of the proposed method and to identify the contribution to the error of the steps of the approach, we have decided to extract predicting equations directly from the measurements, thus evacuating errors related to the modelling steps. Alike Eqs. (10) and (13) that were determined from simulated data, Eqs. (17) and (20) were extracted from the measurements results provided in Fig. 17.
The obtained variational laws for the sensing parameters extracted from measurements are:    (17) and (20), solving the inverse problem from the measured data lead to the values proposed in Table 7.
This table shows a notable improvement in the Relative Difference values for sucrose concentrations and slight improvement concerning sodium chloride Extracted Concentrations with respect to Real Concentrations. It is worthwhile to note that the higher errors presented in Table 7 are related to the same samples than those exhibiting highest error level in Table 6. Improvements in the Relative Differences between Tables 6 and 7 are related to the use of the measured results to extract the equations (Table 7) instead of equations based on the modelling stages (Table 6). This means that the differences between errors presented in Table 7 and those exhibited in Table 6 are related to the modelling stages: solutions characterization and their mathematical modelling on one hand and sensor numerical simulations using these material models on the other hand.
Moreover, the resulting errors presented in Table 7 could be implied by either experimental inaccuracies or related to the regression coefficients during the equation extraction step.
To evaluate the contribution of experimental inaccuracies upon the resulting error, several measurements of some aqueous solution samples were made in a row, all external parameters being kept identical. The resulting relative errors between values extracted from two consecutive measurements appeared to reach up to 3.9%, yielding quadratic differences up to 5.4%. These results corroborate the hypothesis of experimental electrical measurement inaccuracies, which could be attributed to, amongst others, the measurement equipment or the rinsing procedures.
Additionally, errors could also come from the samples themselves, their exact composition and concentrations being unknown. Indeed, as described in "Dielectric characterization and modelling of aqueous solutions" section, samples were prepared by dissolving the ad-hoc weighted compounds under powder form into DI-water. This process might introduce errors, which would be consistent with the fact that errors appear to be consistently on the same samples.

Conclusions
Microwave characterization and modelling of binary solutions involving sodium chloride and water on one hand and sucrose and water on the other hand have been presented, occurring physical phenomena being described. Then, ternary mixtures made of sucrose, sodium chloride and water at various concentrations have (17) |S 11 (100 MHz)| = |S 11 (100 MHz)| H2O − a i C 2 NaCl + b i C NaCl , with a i = 2.29 × C C 12 H 22 O 11 − 2.94, (19) and b i = −0.93 × C C 12 H 22 O 11 + 1.65, (20) and SR Samples = SR DIW − −3.71 × C NaCl × C C 12 H 22 O 11 + 3.65 × C C 12 H 22 O 11 + 5.08 × C NaCl , Table 7. Summary of extracted concentrations using equations derived from measured results (Eqs. (17) and (20)). www.nature.com/scientificreports/ been electromagnetically characterized over a wide frequency band, yielding an EM model parametrized by the sucrose and sodium chloride concentrations. Based upon these models, a bi-parameter microwave sensor able to track both sucrose and sodium chloride in aqueous solutions has been introduced. Its design relies on two electrically-coupled open-loop (SRR) resonators, the sample being located in the inter-coupling area. The sensor design parameters have been optimized during a retrofitting stage to finely adjust the 3D model of the sensor with measurements. This calibration stage was made using air and DI-water as references. Afterwards, sensing parameters have been discussed and chosen, they consist in the reflected S-parameter and a Sensing Ratio relying on the central frequency, bandwidth and transmitted S-parameter, which is calculated around 7 GHz. Then, full-wave simulations of the sensor were run with the previously-described models for samples exhibiting various concentrations of sucrose and sodium chloride. From the simulation results, predictive variation laws of the sensing parameters with respect to the compounds concentrations have been set. Finally, real samples have been measured and the predicted concentrations values have been compared to the real ones, demonstrating the ability of the sensor to retrieve concentrations with a rather good accuracy within the measured ranges.
Further work could include enhancement of the sensitivity through the improvement of the sensor design to address lower concentration ranges, which could be more suitable for potential industrial applications. Finally, the design of a microwave sensor able to monitor more than two compound concentrations is also a worthwhile challenge by considering the craving demand of highly sensitive, continuously-measuring sensors.